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We study T 3 Gowdy spacetimes in the Einstein-Maxwell-dilaton-axion system and show by the 
Fuchsian algorithm that they have in general asymptotically velocity-term dominated singularities. 
The families of the corresponding solutions depend on the maximum number of arbitrary functions. 
Although coupling of the dilaton field with the Maxwell and/or the axion fields corresponds to the 
"potential" which appears in the Hamiltonian of vacuum Bianchi IX spacetimes, our result means 
that the spacetimes do not become Mixmaster necessarily. 



O ' I. INTRODUCTION 

o : 

The singularity theorem gives us some sufficient conditions for the existence of spacetime singularities (a spacetime 



00 



with incomplete geodesies) 19 18|]. These conditions fall into four categories : (1) the strong energy condition; (2) 
Oh 1 the generic condition; (3) the chronology condition and (4) the existence of trapped regions. In addition, after the 
discovery of the singularity theorem, various generalization by replacing these conditions have been suggested by some 
authors (see e.g. Refs. [ pi|p^j3l[ | ) . So, although, we have new many variations of the singularity theorem, none of 
ON \ them gives information on the nature of the singularity. 

To the best of our knowledge, the first research about the nature of the singularity was done by Belinskii, Khalatnikov 
and Lifshitz (BKL) They investigated how the spacetime evolves into the singularity and conjectured that the 
. dynamics of nearby observers would decouple near singularities, (the BKL conjecture.) BKL described a vacuum and 
spatially homogeneous (mainly Bianchi IX) spacetime singularity as an infinite sequence of Kasner epochs ( "oscillatory 
approach to the singularity"). Independently, Misner showed the same behavior in terms of exponentially growing 
"potential" terms constructed by the spatial curvature. The spacetime is bounced by the potentials a infinite number of 
times. Such behavior is called the Mixmaster dynamics (^8|. Then, BKL further speculated that a generic singularity 
should be locally behaved like a Mixmaster type. To date, the BKL conjecture has neither been validated nor 
invalidated in general by rigorous arguments although Ringstrom has recently have a result which supports the 
validity of the BKL conjecture, that is, he has shown rigorouslythat the Bianchi type IX solutions converge to an 
attractor consisting of the closure of the vacuum type II orbits [g2| . 
Ch . There is a special case of the BKL conjecture called the asymptotically velocity-term dominated (AVTD) singularity. 
It is not described by a infinite sequence of Kasner epochs but by only one epoch (i.e., "no oscillatory approach to the 
\ singularity" ) |2~l]j3"6|| . In brief, the characteristic feature of the AVTD singular behavior is that all spatial derivative 
terms of the field equations are negligible sufficiently close to the singularity. Hence, Kasner spacetimes are necessarily 
(A)VTD. It is possible to make rigorous arguments on the nature of the AVTD singularity since the AVTD solutions 
are simpler than the Mixmaster solutions in the sense that there is no complicated oscillation. Our interest is to verify 
BKL conjecture and clarify whether spacetimes is AVTD or not in non- vacuum and/or spatially inhomogeneous cases. 

In order to attack the above issue, we will consider Gowdy spacetimes Jl6| . They are spatially compact spacetimes 
which have U{1) x U(l) symmetry and vanishing twist. Gowdy spacetimes are adopted in various investigations because 
they are one of the most manageable inhomogeneous spacetimes and can give essential features of inhomogeneity in 
spite of its simpleness. The behavior near the singularity in Gowdy spacetime have also been discussed in various 
situations. In the vacuum case, Isenberg and Moncricf have shown that polarized Gowdy spacetimes are AVTD 
Grubisic and Moncrief and Kichenassamy and Rendall |^3| have shown that non-polarized T 3 Gowdy spacetimes 
have AVTD singularities in generic in the sense that the AVTD singular solutions to Einstein's equations depend on the 
maximal number of arbitrary functions if the low velocity condition is satisfied. When one of these functions is constant 
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(i.e. non-generic), it was also shown without the low velocity condition that non-polarized T 3 Gowdy spacetimes are 
AVTD |2^] . Using numerical method, Berger, Garfinkle and Moncrief found that non-polarized T 3 Gowdy spacetimes 
have AVTD singularities even in the situation that the low velocity condition is broken initially Recently, 
Garfinkle has shown that S 2 x S 1 Gowdy spacetimes have AVTD singularities numerically [pH[ . 

AVTD solutions with matter fields present have been found. Such examples include polarized Gowdy spacetimes 
admitting scalar fields minimally coupled with the Maxwell field Jll| and the Einstein-dilaton-axion (EDA) system in 
polarized Gowdy spacetimes It has been shown numerically, however, that magnetic Gowdy spacetimes are not 
AVTD but Mixmaster 0,0. Hence, it is nontrivial to determine whether the spacetime is Mixmaster or AVTD type. 
Because of these facts, we shall study the non- vacuum and spatially inhomogeneous case. 

There is an another point which needs to be clarified. Belinskii and Khalatnikov have suggested that the existence 
of massless scalar fields suppresses Mixmaster behavior by examining the Bianchi I spacetimes It is due to the 
fact that the scalar field is algebraically equivalent to the stiff matter in some cases. Recently, Berger found some 
confirmation for this suggestion by numerical analysis g]. In both papers, however, it was suggested that exponential 
coupling of the scalar field to the Maxwell field restores Mixmaster behavior which was suppressed initially by the 
scalar field dJ^]. It is important to note that such matter fields arises naturally from low energy effective superstring 
theory |35|, i.e. the Einstein-Maxwell-dilaton-axion (EMDA) system. The EMDA system has been discussed actively 
in the context of the black hole solutions and singularities, and gives characteristic features different from the Einstein- 
Maxwell (EM) system @. 

Our purpose in this paper is to explore the nature of singularities in the EMDA system with Gowdy spacetimes on 
T 3 x R. We make use of the Fuchsian algorithm developed by Kichenassamy and Rendall 23 1 . Although the system 
treated here is very complicated, the Fuchsian algorithm is powerful enough for the system since this algorithm is 
independent of non-linearity, number of functions and dimension of space. We will show that there exist families 
of AVTD singular solutions of the EMDA field equations, which are generic in the sense that they depend on the 
maximal number of arbitrary functions. Hence, our result shows that exponential coupling of the scalar field to the 
Maxwell field does not necessarily lead the Mixmaster behavior. 

Organization of this paper is as follows. In Sec. [n[ an appropriate definition of AVTD is given. In Sec. 



Ill 



discuss the importance of studying the EMDA system. Sec. |V] is devoted to introducing T 3 Gowdy spacetimes in 
the EMDA s yste m. Sec. presents a brief review of the Fuchsian algorithm. Our main results are given in Sec. VI. 
Finally, Sec. VII is a summary. 



II. ASYMPTOTICALLY VELOCITY-TERM DOMINATED (AVTD) SINGULARITY 



In this section, we give precise definition of AVTD. In contrast to the Mixmaster singularity which is complicated, 
the AVTD singularity is called simple or quiescent An essential difference between Mixmaster and AVTD 

singularities is whether spatial curvature terms in Einstein's equations becomes dominant or not as the system 
approaches the singularities. As mentioned in the previous section, the "potential" terms constructed by the spatial 
curvatures grow exponentially in the Mixmaster case. These "potential" terms give rise to a infinite number of 
reflection and the spacetime will approach an oscillating state, i.e. Mixmaster singularities described as an infinite 
sequence of Kasner epochs. On the other hand, such a "potential" term does not exist in the AVTD case. Therefore, 
reflection by the potential does not occur and the spacetime will approach a single Kasner epoch. 

Now, let us define AVTD. Suppose that a four-dimensional spacetime (M,g) with the signature ( — h ++) satisfies 

Einstein's equations, G^ v — T^, where G p „ = ^R^v — k 1 

momentum tensor. For any 3 + 1 decomposition, we can obtain the constraint equations, 



Rg^v is the Einstein tensor and T M „ is the energy- 



(3) R _ K ab K ab + {trK f = 2T °, 



(2.1) 



and the evolution equations, 
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where h ab , K ab and (trK) are the first and second fundamental form and mean curvature of the 3-dimensional spacelike 
hypersurface, respectively. Also, ( 3 V, ^R a b and ^R are the spatial covariant derivative, spatial Ricci and scalar 
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curvature, respectively. N and N are the lapse function and shift vector. £n is the Lie derivative along the vector 
field N. From the full Einstein equations, we define the velocity term dominated (VTD) equations as follows: 



K ab K ab + {trKY = 2T U , 



(2.5) 



^ a KS- {3) V b (trK) = -f b °, 



(2.6) 



d t h ab = -2NK ab 



(2.7) 



3 t Kl = N [(trK)K%] + 2N 



f b a + -h a b (f°~f c c 



(2. 



where • is modification of •. (Generally, the spatial derivatives are removed.) 

In this setting we can give a definition of an AVTD spacetime and its singularity. 



Definition 1 (Isenberg and Moncrief [21]) A spacetime {M,g) is AVTD if 

1. it is a solution to Einstein's equations; 

2. there exists another (M, g) and there exists a foliation i t : X 3 — ► M such that 

(a) (M,g) obeys the VTD equations relative to it; 

(b) metric g asymptotically approaches g as t — > in the sense that, for an appropriate norm || • || on the space 
o/3 + 1 quantities •, for any e > there exists T > such that || • — i|| < e for all t < T . 

Furthermore, a spacetime (M,g) has the AVTD singularity if 

1. it contains a spacelike singularity at t — > 0; 

2. it is AVTD. 

To sum up, a singularity is called AVTD if all spatial curvature and spatial derivative terms in Einstein's equations 
can be neglected near the singularity t — > 0. Then, near the AVTD singularity, AVTD solutions can be interpreted as 
a different spatially homogeneous cosmology at each point in space. 

In order to see an example of the AVTD singularity, let us consider vacuum Bianchi I (Kasner) spacetimes. Since 
these spacetimes are spatially homogeneous, there is no spatial derivative terms. Furthermore, they are spatially flat, 
i.e. = 0. Therefore, the singularity is necessarily AVTD. On the other hand, Bianchi IX spacetimes are not 

AVTD but Mixmaster [pp§|] since the spatial curvature does not vanish and the "potential" term constructed by the 
spatial curvature cannot be neglected as approaching to the singularity (See Sec. III). 



III. INFLUENCE OF MATTER FIELDS 



BKL assumed the vacuum spacetimes originally in the investigation of the singularity. The reason for putting such 
assumption is illustrated as follows. Consider class A Bianchi spacetimes with the metric 



ds 2 



-dt 2 



S 2n (e 2 ^) i 



G G J 



where e° is the averaged scale factor, (3 = diag(— 2/3 + , /3 + + V3/3_, /?+ — V3/3-) is a traceless matrix that determines 
the anisotropy and da 1 = e % ik^ A a . For simplicity, take a % = dx l , i.e. Bianchi I spacetimes. Then, e 



3f2 
f-2 



t and 

-60, 



il 2 = fi\ J r0 2 _. The spacetimes have singularities at — > — oo. The square of shear a is estimated as a 2 ~ t~ " ~ e 
Further, consider a perfect fluid whose equation of state is p = (7— l)p, where p and p are the energy density and the 
pressure, respectively. The contribution from the shear is more dominant than that of the energy density p ~ e~ 3 ^ 7 
near singularities (CI — > —00) if 1 < 7 < 2 for natural matters such as a dust (7 = 1) and a radiation (7 = 4/3). 
Thus, the effect of the matter field is negligible for the behavior of the spacetime near singularities, and we have only 
to consider the vacuum case. Indeed, the known Bianchi perfect fluid solutions have similar behavior as the vacuum 
solutions near the singularities Bq]. 
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We must, however, look more carefully into the influence of the energy density in the case 7 = 2, where time 
dependence of the energy density and shear are contribute equally. The matter with 7 = 2 is a stiff matter ||, and 
is not realistic matter in the relativistic sense. However, it is known that a scalar field <f> is algebraically equivalent 
to a stiff matter if V0 is timelike where V is covariant derivative with respect to g. Therefore, the existence of 
scalar fields will give strong influence on the nature of singularities. 

Actually, Belinskii and Khalatnikov found that a massless scalar field can suppress Mixmaster oscillations § . The 
result obtained by Belinskii and Khalatnikov is that all of Kasner exponents can be positive in Bianchi I spacetimes if 
the scalar field exists. Therefore, complicated permutation in the direction of anisotropic contraction toward singular- 
ities does not appear. Numerically, Berger confirmed Belinskii and Khalatnikov's result in Bianchi IX spacetimes and 
furthermore, she showed that Mixmaster oscillations are suppressed in non-polarized J7(l)-symmetric and magnetic 
Gowdy spacetimes which are Mixmaster if there exist no scalar fields (|]. 

On the contrary, there is a possibility that the existence of exponential potentials of scalar fields drastically changes 
the behavior of singularities, i.e. the Mixmaster dynamics may come back. Recall the mechanism of the Mixmaster 
dynamics. The Hamiltonian Hix for Bianchi IX spacetimes ||^|| is 

m IX = -vl+vl+P 2 . + v(n,0±), 

where and p± are conjugate momenta of Q and [3±, respectively. The "potential" term V(Q,/3±) is given by 
V(Q,f3±) = e 4n [e~ 8 ?+ + e ^+W5/J_ + e ^+-^0- _ 2 (e 4/3 + + e ~^~^- + 

In the Bianchi I spacetimes there is no such "potential" term. It is known that Bianchi IX spacetimes have a singularity 
as il — > —00. If the spacetimes are AVTD, the "potential" term must vanish as £1 — > —00. When V(f2, /?±) = 0, we 
have f}± = p± + v±\£l\, pn =const. and p± =const. by varying the Hamiltonian Ttix, where v± = p±/\pn\ and 0± 
are constant. The Hamiltonian constraint TLjx — is in. + v_ =1. Then, 

y ^ e -4|Sl|(l+2cos0) _|_ e -4|n|(l-cos0-V3sin0) _|_ g-4|n| (1-cos 6>+V3 sin 6) 

where we parameterize v+ — cosd and V- — sin 9. For generic 9 (except for 9 = 0, 2tt/3, 4ir/3), the term V expo- 
nentially grows as Q — > — oo. This is a contradiction to vanish the "potential" term as — > —oo. Thus, Bianchi IX 
spacetimes are not AVTD since the "potential" term becomes dominant near the singularity. As a result, reflections 
by the potential cause the Mixmaster dynamics. 

As we have seen in the above discussion, the essential point is the existence of exponentially growing "potential" 
terms. Then, Belinskii and Khalatnikov have claimed that a scalar field exponentially coupled with the Maxwell field 
restores Mixmaster oscillations in homogeneous spacetimes || . Berger also suggested that if there exists a "potential" 

V{4>)=A 2 e c " t ' + B 2 e- l} ' l} , a,P>0, (3.1) 

where A and B are functions, then Mixmaster oscillations can be retained in both homogeneous and inhomogencous 
spacetimes M. 



We wish to explore how "potential" terms such as those in Eq. (3.1) produced by matter fields might change the 
nature of the approach to the singularity. As a first example we consider matter coupled to the Gowdy spacetimc. 
It has been shown mathematically and numerically that vacuum Gowdy spacetimes are AVTD 0,0 23 . It is known 



numerically that certain types of magnetic fields can ruin this behavior |37| . While heuristically one would not expect 
this to happen for the matter fields we study here, no rigorous demonstration that Gowdy coupled to matter fields 
including exponential potentials remains AVTD has been possible until now. 

It should be noted that this form of potential arises from the low energy effective superstring theory. It contains the 
dilaton which is non-minimally coupled with the Maxwell and axion field. From the unified theoretical point of view, 
it is believed that the distinction between gravity and matter fields cannot be maintained near spacetime singularities 
and the superstring theory is the most promising candidate of such unified theories Pq] . Also, it is pointed out that 
the EMD system with a positive cosmological constant has new mechanism forming singularities ^(|. Hence, it is 
important to investigate the behavior around the singularity in the EMDA system. 

IV. EINSTEIN-MAXWELL-DILATON- AXION SYSTEM IN GOWDY SPACETIMES ON T 3 x R 

In this section, we introduce a standard model with additional assumptions. First, we adopt Gowdy spacetimes. 
Gowdy proposed spatially compact and inhomogeneous spacetimes with the following metric M ; 
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ds 2 



e A (t ,0)/2 t - 1/2( _ df 2 + + R ^ e) L-z(t,e) {dy + x ^ e)dz) 2 + e z(t,e) dz s 



(4.1) 



Gowdy spacetimes have two twist free spacelike Killing vectors d/dy and d/dz. The spatial topology of the spacetimes 
is classified into three types, T 3 , S 2 x S 1 and S 3 . We will consider only the simplest case, T 3 , since the isometry group 
action has no degenerate orbits and therefore Einstein's equations have no corresponding singularities p^j. Gowdy 
spacetimes are called polarized if X = 0, and non-polarized if this condition is not meet. 

Properties of the metric (4.1) depend on whether Vi? is timelike, spacelike or null. When the metric Eq. (4.1) 
describes a cosmological model, i.e. Vi? is globally timelike and the spatial topology is T 3 (periodic in 8), one can 
take the function R(t,9) — t without loss of generality by Gowdy's corner theorem if the spacetime is vacuum [ pi] . 
This fact be seen from Einstein's equations, 



Gt 



G ee = R- R" = 0, 



(4.2) 



where dot an d p rime denote t and 9 derivatives, respectively. Also in the EDA system, Eq. (4.2) holds In the EM 
system, Eq. (4.2) is not satisfied generically. However in the case that the field strength F^ v = d^A v — d v A^ has only 
the following components p7fl; 



F, 



ty - U, 



F, 



X, Fez 



X 



Eq. (4.2) is satisfied. This situation does not change even if the dilaton field is exponentially coupled with the Maxwell 
field. Thus, in the EMDA system, we can put i ?(t, 8) = t. Hereafter, we will choose this gauge. In this case, Gowdy 
spacetimes have spacelike singularities at t = |16| . 
The action of the EMDA theory is 



S = 



d A x^ 



1 

--t 
2 



A — A" + Z 2 

z v 2 



WR + e 

- Z 



-2e- 2act > [(e~*X 



-2t 



(0 2 -0' 2 ) 



1 



-2ac 



F 2 + 2 (V0) z 



3 



,2 , e -2Z x2 



e z ) (u 2 



,'2 



X' 



)+e- Z (x 2 -x' 2 )-2A^(c0x-cY)] 



(4.3) 



where <f) = <j>(t, 9) is the dilaton field, H = dB = — ^e 4 "^ * dn(t, 8) is the three-index antisymmetric tensor field dual 
to the axion field n, and a is a coupling constant |3q] . The functions are periodic in 8 (0 < 9 < 2n) because the 
spatial topology is T 3 . [However, the spatial compactness is not relevant to our analysis.] In the case of the effective 
superstring theory, a = 1. If we neglect the axion field, a — \/3 gives Kaluza-Klein theory. If a = 0, the dilaton and 
axion fields are decoupled with the gravity. 

Note that the metric function A is decoupled with other functions, and that A appears only in the Hamiltonian 
and momentum constraints. This is the primary advantage of Gowdy spacetimes. Now, let us focus on evolution 
equations. We can calculate the metric function A by evaluating the integral of A' from to 2ir after obtaining other 
functions. 



V. REVIEW OF THE FUCHSIAN ALGORITHM 



The above system is very complicated and gives rise to highly nonlinear partial differential equations which are 
nontrivially coupled with the matter fields. We would like to construct AVTD singular solutions of such complicated 
equations which are parametrized by as many arbitrary functions as possible. In order to do this, we adopt the Fuchsian 
algorithm developed by Kichenassamy and Rendall j2^| . One of the advantages of this algorithm is applicable to general 
(nonlinear and/or singular) partial differential equations (PDEs). Another advantage of the Fuchsian algorithm is 
applicable to arbitrary spatial dimension. [Spatial dimension in Gowdy spacetimes is one since these spacetimes have 
two spacelike Killing vector fields.] This fact suggests that the Fuchsian algorithm is applicable to spacetimes with 
fewer or no symmetry p^Jl[]. 

In this section, we briefly review the Fuchsian algorithm p|J2(lp3]. Let us consider a PDE system 

F[u{t,x a )} = 0. (5.1) 

Generically, u can have any number of components. Here, we will assume that the PDE is singular with respect to 
the argument t. The Fuchsian algorithm consists of three steps as follows: 
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Step 1 Identify the leading (singular) terms uo(t, x) which are parts of the de sire d expansion for u. This means that 



the most singular terms cancel each other when u${t, x) is substituted in Eq. (5.1). 
Step 2 Introduce a renormalized unknown function v(t, x), which is given by 

u = u Q + t m v. (5.2) 

If uq ~ t k , we should set m = k + e, where e > 0. Thus, v is a reg ular par t of the desired expansion for u. 
Step 3 Obtain a Fuchsian system for v by substituting Eq. ( |5.2| ) in Eq. (5.1). That is, 



(D + A)v = t e f(t,x,v,d x v), (5.3) 
where D = i<9( and A is a matrix which is independent of t and e > 0. / can be assumed to be analytic in all of 



arguments except t and continuous in t. Note that Eq. (5.3) is a (singular) PDE system for the regular function v. 



Roughly speaking, the Fuchsian algorithm is a transformation fro m fin ding singular solutions of the original equa- 



tions (5.1) into finding regular solutions of the Fuchsian equations (5J3) which may be singular even if the original 
equations are regular. Once we have the Fuchsian system, we can show the existence of a unique solution for prescribed 
singular part uq by the following theorem. 



Theorem 1 (Theorem 3 of Ref. [23]) Let us consider a system Eq. (5.S), where A is an analytic matrix near 



x = 0, such that \\o~ || < C f or < a < 1 (boundedness condition) and f is analytic in space x and continuous in 



time t. Then the system (5.1 ) has exactly one solution which is defined near x — and t = 0, and which is analytic 



in space and continuous in time, and tend to zero as t — ► 0. □ 

Applying the Fuchsian algorithm to our problem of Gowdy spacetimes in the EMDA system, the procedure is 
summarized as follows: 

1. Solve VTD equations as ordinary differential equations with respect to t. Solutions obtained such a way are the 
leading (singular) terms of the desired formal solutions. 

2. Substitute the formal solutions into full field equations. Here, make integral constants arbitrary functions of 
spatial arguments x. 

3. Obtain a Fuchsian system for unknown functions and evaluate every eigenvalues of A. 

4. Apply Theorem |j] to the system. 

If every eigenvalues of A are non-negative, the boundedness condition in Theorem [j] holds. Then, renormalized 
unknown functions must vanish as t — ► by Theorem [jj. Therefore, the only singular terms which are solutions to 
VTD equations remain and they are solutions to full field equations at t = 0. Thus, we have AVTD singular solutions. 
Essentially, the above argument is a singular version of the Cauchy-Kowalevskaya theorem Jl|,^3| . 

VI. ASYMPTOTIC BEHAVIOR IN GOWDY SPACETIMES 

We are now ready to examine the singularity in the EMDA system. The EMDA system with Gowdy symmetry, 
however, has six functions (Z ', X, u>, \, 4>, ^), and it is very complicated to treat all of them at the same time. Hence, 
we will reduce the system by selecting four functions from them for simplicity. As a selection rule, it is important to 



leave exponential coupling between the dilaton and other fields since, as mentioned in Sec. Ill, an essential part of our 
analysis is whether or not "potential" terms in Einstein's equations exponentially grow as t Therefore, we 

will consider three cases, which are the non-polarized EMD, the non-polarized EDA and the polarized EMDA. 

A. non-polarized EMD (x = 0, k = 0) 

First, let us consider the case that the dilaton field is exponentially coupled with the Maxwell field. Varying the 



action (4.3) with respect to (Z, X, cj, x), we obtain evolution equations, 



G 



D Z Z - t z Z" + t 2 e- 2Z (X 2 - X' 2 ) + 2te- 2a *(e- z X' 2 - e z )(dj 2 - J 2 ) = 0, 
D 2 X - t 2 X" - 2t 2 (XZ - X'Z') - Ate- 2a ' l ' +z X(6j 2 - J 2 ) = 0, 



D 2 uj - t 2 cu" - t 2 



+t 2 



■ 1 {2X - ZX)Xer z + Ze z 
2 «^+7 e -z x2 + e z 

<2X' - Z'X)Xe- z + Z'e z 



D 



t 



t 

2acj)' 

2 ill 



ate 



e- z X 2 

-Z v 2 



J = 0, 



(e- z X 2 + e z ) {uo 2 -J 2 )=Q. 



(6.1) 
(6.2) 



(6.3) 
(6.4) 



Following the pr ocedure e xplained in the previous section, we have VTD equations by dropping spatial derivative 
terms from Eqs. (|6.l| )-( |6.4| ). 



D 2 Z + t 2 e~ 2Z X 2 



D^uj - r i 2acf) + - - 



2te- 2a,p (e- z X 2 - e z )uj 2 = 0, 
D 2 X - 2t 2 XZ - Ater 211 ^ 2 Xlo 2 = 0, 
1 (2X - ZX)Xer z + Ze z ^ 



'X 2 



•w = 0, 



+ e' 

-Z v 2 



(6.5) 
(6.6) 

(6.7) 
(6.8) 

According to Definition [j], solutions to these equations are AVTD if they exist and are also solutions to full Einstein's 
equations as t — > 0. 

Since we expect Kasner-likc solutions, Z w Zq hit, X w Xq and <p <f)Q lni are chosen as the leading terms for Z , X 
and cj>, respectively. Substituting them into Eq. (6.7), we have a leading term of u). Thus, we can construct a family 
of formal solutions to full Einstein's equations as, 



D 2 (j) + ate- 2a ' ) '(e- z X 2 + e z )uj 2 = 0. 



Z = Z (d) \nt + Z x {9) + t e a(t, 9), 
X = X o (0) + t 2Zo {X 1 (9) + 0{t,6)) , 
w = loq{9) + t z "+ 2a ^' +1 ( Wl ((9) + 7 (i, 9)) , 
= 0o(6 ) ) lnt + <M0) +t e 6(t, 9), 



(6.9) 
(6.10) 
(6.11) 
(6.12) 



where Zq > 0, Z + 2a0 o + 1 > 0, e > 0. The first and second inequalities are given by Eqs. ( |6.6|) a nd (3.7j ). a, j3, 
7 and 5 are renormalizcd unknown functions. Note that when we put a — ( 3 = 7 = 5 = 0, Eqs. (5.9)-(3.12) are exact 
s olu ti ons to VTD equations as t — > 0. Substituting the formal solutions (6.S)-( 6.12 ) into full Einstein's equations 
( 3.1 )-( 3.4 ), we can obtain the following system, 



(D + A)u = f(t,9,u,it), 

where 

ii= (a, Da, ta', /3, D/3, 0, 7, Dj, tf, 5, D5, tS') T . 
The matrix A has the 12 x 12 components, 



A = 



(6.13) 
(6.14) 




where, 



A z = 





-1 






(0 


-1 




{ 


2e 


• 


A x = 





2Z 


• 





















z 





-1 

2a<po 






1 




As = 



0-10 
e 2 2e 
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We can easily verify that eigenvalues of A are 0, 2Zq, Zq + 2acj)Q + 1 and e, which are non- negative. Then, we conclude 
that the boundedness condition of Theorem |l| holds. Next, we examine the leading parts in components of the vector 

/ \ 

X l2 e -2(Z 1 +t'a) t 2-2Z Q -t _|_ 2X 2 UJq e — ( 2a 1 +*" &) + Zi+t' a) ^1- Z o -2a0 o -e _| 

t(a + Da)' 


-2X' (Z' lnt + Z[ + t e a)t 2 - 2Zo - ^^e-^'^+^'+^+'^i 1 ^ -^ + ■ ■ ■ 

tip+Dpy 



X- 2 e Zl + t€a {uj% - [2a(4>' lnt + cj)[ + t e 8') + lK}* 1 "^- 2 ^ + • • • 



a ^2 u; /2 e -(2a(0i+t e (5) + Z 1 +t e a) i l-Zo-2a0o-e _| 

V " 1 t(8+Dsy ) 

The dots represent the terms which are lower order with respect to t and are not important for our discussion since 
the power of these terms is necessarily positive if conditions Zq > and Zq + 2a<j>o + 1 > hold. Positivity of the 
power of t in components of the vector / gives us sufficient conditions to obtain the AVTD solutions. Finally, we 
apply Theorem [l] to this system and then, we obtain the following theorem. 

Theorem 2 Let Zi(ff), Xi(9), tUi(0), 4>i(9), where i = 0, 1 , be r eal analytic functions and < Zq < 1 and — 1 < 
Zq + 2acf)Q < 1 for <0< 2tt. Then, for field equations \6. in the EMD system, there exists a unique solution 
with the form \6.S[ )-( 6A^ ), where a, [3, 7 and 8 tend to zero as t — > 0. □ 

Theorem || means that T 3 Gowdy spacetimes in the EMD system have AVTD singularities in "general" in the sense 
that the singular solutions have maximal number (i.e. eight) of arbitrary functions (Zq, Xq, lvq, <j)Q, Z\, X\, w\, (j>\). 

In components of the vector /, terms t 1 - z o- 2a <t>o anc l ^2-2Z gj ve us U pp 6 r bound for Zq and Zo + 2a<po in Theorem||. 
However, since these terms are always multiplied by Xq and lo'q, the upper bound can be removed if Xq = and 
lu q = 0. Hence, we obtain the following corollary: 

Corollary 1 Let Zi(9), Xi{9), uJi{6), 4>i(9), where i — 0, 1, be real a naly ti c fu nctions, and assume < Zq, —1 < 
Zq + 2a<f)Q and Xq = uj' = for < 9 < 2ir. Then, for field equations (6.1 )- j\6-4 ) in the EMD system, there exists a 
unique solution with the form ^6.!\ )-(6.1i ), where a, (3, 7 and 8 tend to zero as t — > 0. □ 

Note that this is a "non-generic" case in the sense that two of arbitrary eight functions are constant in 9. 



B. non-polarized EDA (uj = 0, \ — 0) 



Now, we will consider the case that the axion field is exponentiall y c oupled with the dilaton field. In the similar 
way to Sec. VIA, we have evolution equations by varying the action (4.3) with respect to (Z, X, <f>, k). 



D 2 Z - t 2 Z" + t 2 e~ 2Z (X 2 - X' 2 ) = 0, 
D 2 X - t 2 X" - 2t 2 (XZ - X'Z') = 0, 
D 2 cf> - t 2 4>" - -t 2 e^(k 2 - k 12 ) = 0, 



D 2 K 



t 2 n" + 4at 2 (<j>k- 4>'k') =0. 



Note that a pair of PDE s for gravi ty ( 6.15 ) and ( 6.16 ) is decouple d wi th tho se for scalar fields ( 6.17 ) and 
Furthermore, Eqs. (6.15) and (6.16) have equivalent form to Eqs. (6.17) and (6.18) as simultaneous PDEs. 
facts are known as "mirror images" in string cosmologies and used to construct new nontrivial solutions [p4| . 
Dropping the spatial derivative terms from Eqs. ( 6.15| )-( |6~T^ ), VTD equations are obtained as follows: 



D 2 Z + t 2 e- 2Z X 2 = 0, 
D 2 X - 2t 2 XZ = 0, 

D 2 - -tV^A 3 = 0, 
2 

D 2 K + 4ai 2 (/>K = 0. 



(6.15) 
(6.16) 
(6.17) 

(6.18) 

(B. 

These 



(6.19) 
(6.20) 
(6.21) 

(6.22) 



8 



We can easily solve these equations as t — > and find formal solutions to full Einstein's equations. 



Z = Zq (9) In t + Zi (9) + t e a(t , 9) , 
X = X (9) + t 2Z °(X 1 (9) + p(t,6)), 

( j ) = MQ)^t + M0) + t e j(t,9) 1 

k = k„(0) + t- 4a ^ ( Kl (0) + 5(t, 9)), 



(6.23) 
(6.24) 
(6.25) 
(6.26) 



where Zq > 0, — 2a0o > and e > 0. These restrictions for Zq and 4>q given by Eqs. ( 6.2C ) and ( 3.22), respectively. 

The system for regular functions a, [3, 7 and 5 in the EDA system is given by Eqs. ( 6.13| ) and (3.14) by replacing 
the 12 x 12 matrix A by 



A = 




where 




A K = 



Eigenvalues of A are 0, 2Zq, — 4a0o and e, and all of them are non-negative. These facts imply that the boundedness 
condition in Theorem [l] holds. After long calculation, the leading parts in components of the vector / is obtained as 
follows. 





x a e -2(z 1+ t'- a ) t 



2-2Z -e 



\ 



t(a + Da)' 


[Xg - 2X (Z' In* + Z' x + t*a')\ t 2 - 2Z ° + ■ 
t(/3 + D0)' 


_S K ' 2 g4a(0i+t' ! 7)^2+4a0o-e -)-... 

t<n + Dj)' 


t(s + Dsy 



J 



Since the power of t of other terms in / is positive under conditions Zq > and —2a(f>o > 0, and then they are not 
needed here. The AVTD behavior requires that the power of the leading parts of t is positive. Thus, we have the 
following theorem which states that non-polarized Gowdy spacetimes in the EDA system are AVTD in general. 

Theorem 3 Let Zi(9), Xi(9), 4>i{9), k, A9), w h ere i — 0, 1 7 be real analytic and < Zq < 1 and < — 2a<f)Q < 1 for 
< 9 < 2tt. Then, for field equations ( 6.1f\) -( 6.1l\ ) in the EDA system, there exists a unique solution of the form 
(6.2c )- (6.26 ), where a, [3, 7 and S tend to zero as t — > 0. 



□ 



As seen from the vector /, the upper bound for Zq and — </>o in Theorem || are derived from terms multiplied by X'q 
and k' q . Thus, we obtain the following corollary in the "non-generic" case. 

Corollary 2 Let Zi(9), Xi(9), 4>i(9), Ki(9), where i = 0, 1, be r e al ana lytic, and assume < Zq, < — 2a0o and 
X' n = k[, = 0, for < 9 < 2ir. Then, for field equations ( 6.1t\ j-(6.1S ) in the EDA system, there exists a unique 

□ 



lution of the form ( 6. 2i\ )- ^6.2t\ j, where a, [3, 7 and 8 tend to zero as t — > 



C. polarized EMDA (X = 0, x = 0) 



Finally, we will study the case of polarized Gowdy spacetimes in the EMDA system. As mentioned in Sec. Ill 



the existence of "potential" like Eq. (|3.l[) might lead to Mixmaster oscillations. Such a potential appears in the 
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EMDA system [See Eq. 3.2£]. On the other hand, it has been shown that vacuum polarized T 3 Gowdy spacetimes 
are AVTD Q. Which behavior is realized in the EMDA system, AVTD or Mixmaster? 

Varying the action (fh3t) with respect to four functions (Z,lj, <f>, k), we have evolution equations as follows: 



d z z -rz 



2 ryll 



2te 



-2a<t>+Z rjji 



0, 



D 2 uj + t 2 



Z - 2a(j) - - J Cj - t 2 u" - t 2 (-2a<j)' + Z')J = 0, 

j-,2 j ,2i« i . -2a<b+Z t ■ 2 I2\ , 1 ,2 4a<ft/-2 I2\ n 

u <p — t tp + are v [u> — u> ) + —at e v (k — k ) = V, 



(6.27) 
(6.28) 



(6.29) 

D 2 k - tV + Aat 2 (kj) - n'cf/) = 0. (6.30) 
Dropping the spatial derivative terms from Eqs. (6.27)-( 6.30| ), following VTD equations are obtained as follows: 

D 2 Z - 2te- 2acl > +z uj 2 = 0, (6.31) 

(6.32) 



D 2 uj + t 2 I Z - 2a<h - - ) Cj = 0, 



D 2 <\> + ate- 2a * +z Co 2 + -at 2 e Aa *k 2 = 0, 
D 2 k + Aat 2 k(j) = 0. 



(6.33) 
(6.34) 



Similarly to the case of the non-polarized EDA system, we can solve these equations as t — > 0. Thus, a family of 
formal solutions of this system can be obtained as follows: 



Z = Z o (0) Xnt + Z^O) +t e a(t,9) 
uj = uj (6) + r 2a ^- Zo+1 ( Wl (0) + 0(t, 6)) 

k = k o (0) + t- ia *° (ki(6) + S(t, 6)) 



(6.35) 
(6.36) 
(6.37) 
(6.38) 



wher e 2a0 p — Zq + 1 > 0, — 4 a^o > and e > 0. Eqs. (6.32) and Q6.34 ) restrict the values of Zq and 4>o ■ By s ubstit uting 
Eqs. ( |6.35| )-( |6.38| ) into Eqs. ( |6.27| )-(|6.30|), the system in the polarized EMDA system is given by Eqs. ( |6.13| ) and fl6.14j ) 
with A given by 



where 







A = 



-1 



( A z 
A u 
A 
V A K 







0-10 



A w = 2a^ - Z + 1 , A K = I -4a0 o 























We can easily evaluate the eigenvalues of the matrix A as 0, 2a0o — Zq + 1, —Aa<j)Q and e. Since they are all non- 
negative, the boundedness condition in Theorem |l| hold. The leading parts in components of the vector / are as 
follows: 



/ = 



/ 

-2ld' 2 ! e -2aWi+*' z l) + Zi+t' 'a t l-2act>o+Z -<i _| 

t(a + Da)' 


H' + u' [-2a(cj)' \nt + <\>\ + r e 7 ') + Z' Q Int + Z[ + t e a']}t 1 - 2a 'l'°+ Zo +■■■ 

t(P + D0)' 


a(jJ l2 e -2a(<p 1 +t'-y)+Z 1 +t' afl-2a<j) +Z -e i a ^12 e Aa{4, 1 +t e )-j-2+ia^> -t _)_... 

r( 7 + ^7)' 2 ° 


[«{,' + 4art' (^ lnt + cj}[ + r £ 7 ')] t 2+4a *» + ■■■ 

t(s + Dsy 
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We do not need other terms since their power is positive under the conditions 2a0o — Zq + 1 > and — 4a0o > 0. 
Thus, we have again the theorem which states t hat p olarized Gowdy spacetimes in the EMDA system are AVTD in 
general even if there exists a "potential" as Eq. (3.1) in the system. This result is a negative answer of the question 



in Sec. III. 



Theorem 4 Let Zi(9), u>i(0), 4>i{9), Ki(9), where i — , 1, be real analytic, and — 1 < 2a0o — Zq < 1 and — 1 < 
2a<po < for < 8 < 2tt. Then, for field equations (6.21)-(6.3L ) in the EMDA system, there exists a unique solution 
of the form ( 6.3b\ -( 6. 3(\ ), where a, (3, 7 and 5 tend to zero as t — > 0. □ 



The leading terms in components of the vector / always include factors lu' and/or k' . Then, we can be weaken the 
conditions for Zq and 4>o as similar to the cases of the EMD and the EDA systems if we put oj' — n' = 0. 

Corollary 3 Let Zi(9), LOi(9), 4>i(9), Ki(9), where i — 0, 1, be re al analy tic, and assume u>' = k' = 0, — 1 < 2acf>o — Zq 



and 2a<f>Q < for < 9 < 2ix . Then, for field equations (6.2 r i)-(6.3C) in the EMDA system, there exists a unique 
solution of the form (6. 3!\)- (6. Sty, where a, {3, 7 and 5 tend to zero as t — > 0. 



□ 



VII. SUMMARY 

Theorem |[ || and [| state that T 3 Gowdy spacetimes in the EMDA system have AVTD singularities in "general" in 
the sense that these AVTD singular solutions depend on the maximal number of singular data. In the "non-generic" 
cases where some functions in the singular data are constant in 9, we have obtained corollary |l], ^ and [3| These 
results support BKL conjecture, that is, the dynamics of spatially different points effectively are decoupled from each 
other. However the spacetimes in our system do not show an oscillation of the Mixmaster spacetime near cosmological 
singularities for some parameter regions. In this sense, our system is a special case. 

When we take the matter fields into account, Theorem^ and corollary || give an answer, "NO", to the question in 



Sec. [II, i.e. the existence of a "potential" Eq. (3.1) does not necessarily restore the Mixmaster behavior. We must, 
however, impose another condition on the initial matter field in addition to the low-velocity condition to realize the 
AVTD behavior by the effect of the "potential". Thus, we can say that the AVTD behavior is controlled by the 
"potential" . 

We can verify that the Kretschemann invariant R^ vXa R^ v \a of all of our AVTD solutions blows up as i - (- z o+ 2 ^o+ 3 ). 
Thus, the AVTD singularity in the EMDA system is the curvature singularity. This supports the validity of the strong 
cosmic censorship for T 3 Gowdy spacetimes in the EMDA system since the spacetimes cannot be extended beyond 
the singularity similar to the vacuum case p?| , ^l| , p9| , ^2[ . 

Let us compare our results with those obtained by Grubisic-Moncrief and Kichenassamy-Rendall in the (non-) 
polarized and vacuum case (T^j2^|. The low- velocity condition derived by Grubisic-Moncrief and Kichenassamy- 
Rendall are equivalent to our condition < Zq < 1. If spacetimes are vacuum, then non-polarized T 3 Gowdy 
spacetimes are AVTD under this condition. Also, when X = (polarized), the spacetimes are AVTD without the 
low-velocity condition. In any case, our results clearly include Grubisic-Moncrief and Kichcnassamy-RendalPs when 
the matter fields vanish. 

We should mention about some recent works related to our analysis. Weaver, Isenberg and Bcrgcr have shown 
numerically that a Gowdy spacetime in the EM system (a magnetic Gowdy spacetime) has the Mixmaster behavior |}7) . 
The magnetic field they considered couples the metric function A through its amplitude. Our Maxwell fields, however, 
do not have such character (cf. [^7|), so Gowdy spacetimes in our EM system are able to show the AVTD behavior. 
It is expected that if Gowdy spacetimes in the EM system such as Weaver's model have Mixmaster behavior, they 
cannot show the AVTD even if we put the dilaton and the axion fields. 

Andersson and Rendall have shown by using the Fuchsian algorithm that a generic cosmological spacetime in the 
Einstein-scalar field system has AVTD property Qj. They have imposed the Gaussian coordinate conditions, goo = — 1 
and go a = on the spacetime. As it is understood from its definition, AVTD is defined under a coordinate condition. 



Then, choice of gauge or coordinate conditions is important. Clearly, our gauge conditions (see Sec IV) are different 
from Andersson-Rendall's. Furthermore, the equations of matter fields used in Ref. Q are linear PDEs. Contrary, 
our field equations in the EMDA system are nonlinear. These nonlinear PDEs cannot be came to linear ones. Thus, 
their and our results can be complements each of the other. 

Damour and Henneaux have shown under the Gaussian coordinate conditions that a generic cosmological singularity 
of low energy effective superstring theory is the Mixmaster [T^ j . Their model is closely related to our model in the 
sense that it contains the matter fields with exponential potentials. However, there are some different points beside the 
gauge conditions : (i) the spacetime dimension (it is 11 or 10 in Ref. Jll|), (ii) an assumption to find formal solutions, 
and (iii) components of p-forms (i.e. Maxwell and axion) fields. It is known that the existence of the Mixmaster 



11 



behavior strongly depends on the spacetime dimension. Also, the contributions of the p-form fields are neglected 
in the field equations when formal solutions, i.e. VTD solutions, are constructed. Then, these formal solutions are 
different from ours since the velocity terms of any fields are not neglected in our analysis. Furthermore, thanks to 
Gowdy symmetry in our analysis, components of p-form fields are restricted, that is, the Maxwell field strength can 



have only four components (see Sec. IV), although there is no such restrictions in Ref. |L3]. 

To avoid complicated calculation we focused only on four functions among the six functions (Z, X, to, <f> and 
k) in the present paper. For a next step, we would like to investigate the full system where that all of six functions 
do not vanish. Another subject to consider is asymptotic behavior of J7(l)-symmetric spacetimes |3(i[| which have 
only one spacelike Killing vector. Although Einstein's equations for [/(^-symmetric spacetimes are more complicated 
than those of T 3 Gowdy spacetime, the Fuchsian algorithm is applicable to such system since this algorithm is 
independent on dimension of space. Recently, the Fuchsian algorithm is used for vacuum polarized t/(l)-symmetric 
spacetimes po|,p|. Our interest is to know whether non-polarized [/(l)-symmetric and generic cosmological spacetimes 
in the EMDA system are AVTD or not. Although the procedure of the Fuchsian algorithm is routine, finding out 
plausible formal solutions to VTD equations and obtaining the Fuchsian system are algebraically complicated, and 
not so easy. Therefore, these subjects are devoted to the future investigation. 
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